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Progress  Summary 

In  the  first  year  of  this  project  (May.  1981  -  April,  1982),  our  efforts 
were  concentrated  on  the  reanalysis  of  the  B-*X  transitions  of  HgCl,  HgBr.  and 
Hgl,  working  with  single  isotopomers  containing  These  transitions  were 
vibrationally  reanalyzed  in  all  three  molecules,  leading  to  significant  changes 
in  the  spectroscopic  parameters,  particularly  for  the  ground  (X  E  )  states. 
At  the  same  time  work  was  begun  on  the  rotational  analyses  of  HgCl  and  HgBr. 
While  these  results  gave  an  Improved  description  of  the  B(  E  )  and  X  states, 
they  were  not  really  adequate  for  the  latter.  The  reason  was  that  the  B-*X 
transitions  do  not  sample  the  low-v  (<7)  regions  of  the  X  states,  because  these 
levels  are  out  of  the  Franck-Condon  region  for  emission  from  the  low  v '  levels 
that  are  mainly  populated  in  our  Tesla  discharge  sources  (and  in  HgX  lasers) . 

To  obtain  spectroscopic  constants  for  low  levels  of  the  X  states,  we 
originally  considered  transient  B*-X  absorption  experiments.  However  these 
experiments  would  have  required  equipment  that  we  did  not  have  and  were  unable 
to  fabricate  in  the  time  available.  Consequently  we  turned  our  attention  to 
the  C(  ^3/2^"*^^  transitions,  which  occur  in  the  UV  for  all  three 
HgX  molecules,  and  which  were  present  In  the  emission  from  our  sources  with 
sufficient  intensity  to  permit  photographing  them  at  high  resolution.  In  every 
case  these  transitions  show  violet-degraded  band  structure  and  terminate  on  low 
V  levels  of  the  X  state.  Therefore,  by  combining  results  for  these  systems 
with  those  for  B-*X,  we  have  been  able  to  obtain  improved  constants,  valid 
typically  for  v  =  0-35  of  the  X  state.  We  have  prepared  a  preliminary  report 
of  this  work  on  HgBr,  which  also  describes  a  new  method  we  have  used  for 
fitting  the  data.  This  paper  is  included  as  an  appendix  to  the  present  report. 
We  have  more  or  less  completed  work  on  the  I>*X  systems  in  HgCl  and  HgBr,  and 
the  C-*X  and  several  weaker  systems  in  Hgl;  and  we  expect  to  prepare  papers  for 
publication  soon.  The  C-*X  systems  in  HgCl  and  HgBr,  and  D-*X  in  Hgl  are  still 
under  investigation.  Work  continues  also  on  the  rotational  analyses  for  HgCl 
and  HgBr. 


I  mentioned  the  development  of  a  new  method  for  representing  vibrational 
energy  levels  in  molecules  like  HgBr(X).  where  the  data  sample  a  large  fraction 
of  the  total  binding  energy  of  the  state.  This  method  is  the  use  of  a  mixed 
representation  for  the  levels  —  the  customary  polynomials  in  (v+1/2)  for  low 
levels,  near-dlssociatlon  expansions  [which  involve  polynomials  in  (▼[)“▼). 
where  v^  is  the  noninteger  vibrational  quantum  number  at  the  dissociation 
limit]  for  high  v.  In  conjunction  with  the  development  of  this  method.  I 
earlier  devised  computational  schemes  for  direct  fitting  of  spectroscopic  data 
to  near-dlssociatlon  expansions,  with  application  to  the  D'->A'  transitions  in 
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Ig  and  BPg  and  the  B-*X  and  D-»X  systems  in  XeCl.  This  work  has  now  been 

published  and  is  included  as  an  appendix  to  this  report. 

Another  new  method  we  have  developed  is  the  use  of  Morse-RKR  potential 

curves.  We  have  found  that  the  repulsive  (left)  branches  of  many  known 

diatomic  potential  curves  can  be  very  well  approximated  by  a  Morse  curve 

derived  from  the  known  R  ,  u  .  and  w  x  values.  For  heavy  diatomics  such  as 

e  e  e  e 

the  HgX  molecules  it  may  be  relatively  easy  to  obtain  vibrational  constants  but 
hard  to  obtain  rotational  constants,  because  the  spectra  are  so  congested.  In 
such  cases  one  can  with  fair  reliability  approximate  the  repulsive  branch  of 
the  potential  by  a  Morse  curve  (guessing  the  unknown  R^) ,  then  obtain  the 
attractive  branch  from  the  RKR  f  integrals  (which  are  calculated  from  the 
vibrational  constants  alone) .  A  paper  describing  this  method  has  been 

published  and  is  included  as  an  appendix  to  this  report. 

Future  Work 

Work  continues  on  the  rotational  analyses  of  B'^X  in  HgBr  and  HgCl.  We 
have  succeeded  in  computerizing  our  line-measuring  procedures  and  are  preparing 
to  apply  our  new  instrumentation  to  the  line-rich  HgCl  spectrum,  and  possibly 
to  HgBr  (where,  however,  line  congestion  makes  necessary  more  oprrator 
intervention).  We  have  also  now  (finally)  set  up  and  begun  testing  our 

Fabry-Perot  interferometer.  We  hope  to  begin  using  it  to  resolve  blended  lines 
and  to  measure  line  widths  as  a  function  of  pressure  within  a  few  weeks.  We 
have  also  recorded  HgBr  and  Hgl  emission  spectra  at  low  resolution  as  a 
function  of  pressure.  We  will  use  these  and  additional  data  to  extract 
extimates  of  the  R-dependence  of  the  B-»X  transition  strength  function. 

Personnel 

Three  graduate  students  —  K.  S.  Viswanathan.  J.  Gail  Ashmore,  and 
0.  Carlysle  Salter  —  have  been  employed  essentially  full  time  on  this  project. 
Dr.  Viswanathan  has  now  completed  the  requirements  for  his  Ph.  D.  degree 
(thesis:  "Part  1:  Spectroscopic  Studies  of  Charge  Transfer  Transitions  in 

Iodine.  Bromine,  and  Mercury  Iodide;  Part  2:  Nitrogen  Laser  Pumped  Dye  Laser." 

June,  1983),  and  has  taken  a  postdoctoral  appointment  in  the  Department  of 
Chemistry  at  the  University  of  Indiana.  Ashmore  and  Salter  will  soon  begin 
their  fourth  years,  but  Salter  will  be  supported  under  another  project.  My 
wife,  Patricia  C.  Telllnghuisen,  was  employed  one-fourth  time  until  March, 


1;.’83.  at  which  time  she  began  to  be  supported  under  another  project  (but  funded 
as  an  extension  of  this  project;  see  below). 


Or  her  Support 

I  have  received  funding  from  DARPA  for  a  project  called  "XeF 

Spectroscopy."  which  is  funded  as  an  extension  to  the  present  project  for  the 

P' riod  March  1.  1983  -  November  30.  1983.  I  have  also  obtained  funding  from 

tie  Air  Force  Office  of  Scientlc  Research  for  "Spectroscopic  Studies  of  the 

llilogens."  for  the  period  April  1,  1983  -  March  31.  1986. 
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The  B—*X  spectra  of  HgQ  and  Hgl  are  studied  at  high  resolution  for  the  single  isotoi»c  species, 
*®Hg”CI,  *"Hg*‘T,  and  “®Hg'”I.  For  Hgl  the  analysis  indicates  that  the  v’  numbering  should 
be  decreased  by  one  imit  from  the  previous  assignment.  For  both  molecules  the  analyses  deviate 
progressively  from  the  previous  assignments  at  high  v",  extrapolating  to  lower  estimates  of  the 
ground-state  dissociation  energies.  Franck-Condon  calculations  yield  AR,  ( =  R,'  —  A, ') 

=  0.60A  for  HgCl  and  0.49  A  for  Hgl.  The  strongest  laser  features  previously  reported  for  HgCl 
occur  near  the  heads  of  the  overlapped  0-22, 1-23, 2-24,  and  3-2S  t^ds.  The  Hgl  laser  operates 
in  the  region  of  the  0-14, 0-lS,  1-15, 1-16, 2-17,  and  2-18  bands. 

PACS  numbers:  33.20.Kf,  33.10.Gx,  33.70.  -  w,  42.55.Hq 


In  the  search  for  new,  high-power,  UV-visible  lasers, 
oonsiderableattentionhasbeenfocusedontheR— transi¬ 
tions  in  the  diatomic  mercury  halides.  To  refine  the  spectros¬ 
copic  characterization  of  these  systems,  we  have  been  study¬ 
ing  the  emission  spectra  of  isotopically  pure  HgX  molecules. 
We  reported  recently  preliminary  results  of  our  work  on 
HgBr.*  In  this  letter  we  discuss  the  results  of  our  vibrational 
analyses  for  HgCl  and  Hgl. 

The  emission  spectra  were  obtained  using  equipment 
and  procedures  similar  to  those  described  previously.  The 
tesla  discharge  sources  were  charged  initiaUy  with  the  de¬ 
sired  isotopic  HgX  species,  which  were  prepared  in  situ  as 
described  for  HgBr. '  Most  of  the  work  involved  sources  con¬ 
taining  ^®“Hg  ”C1  and  ^’‘’Hg  ‘*’1.  To  determine  the  vibra¬ 
tional  numbering  we  also  photographed  and  measured  spec¬ 
tra  of  ““Hg  ( '”l2  from  Oak  Ridge,  stated  isotopic  purity 
99%),  and  *®“H^’C1  (from  a  source  made  with  natural  Clj). 
HgCl  spectra  were  photographed  over  the  region  4400-5800 
A,  at  a  reciprocal  dispersion  of — 5.3  A/mm.  The  Hgl  spec¬ 
tra  were  recorded  for  the  region  4050-4500  A;  most  of  the 
assignments  were  from  plates  having  a  reciprocal  dispersion 
of  ~  1.1  A/mm,  with  a  few  obtained  from  plates  exposed  at 
the  lower  resolution  employed  for  HgO. 

As  in  the  case  of  our  earlier  work  on  HgBr,  our  interpre¬ 
tation  of  the  spectra  of  both  molecules  agrees  qualitatively 
with  WielandV'^  previous  work,  except  at  the  long-wave¬ 
length  end.  However,  our  least-squares  analysis  (see  below) 
indicates  that  the  v’  numbering  suggested  by  Wieland  for 
Hgl  should  be  decreased  by  one  unit.  (Wieland’s  numbering 
was  stated  to  be  uncertain  for  lack  of  a  halogen  isotope  effect, 
since  natural  I2  is  100%  '^^12-)  Below  v’  =  14  in  Hgl  and 
u'  =  20  in  HgQ,  our  measurements  lie  uniformly  2-3  cm~ ' 
below  Wieland’s,  which  is  consistent  with  the  small  isotope 
shift  for  ^‘’"HgX  versus  the  ^HgX  which  dominated  Wie¬ 
land’s  spectra.  However,  for  higher  v’  our  assignments  de¬ 
viate  progressively  to  the  blue  of  Wieland’s,  with  the  discrep¬ 
ancy  amounting  to  30  cm~'  for  our  highest  presently 
assigned  level  (v*  s  23)  in  Hgl,  and  about  40  cm~ '  in  HgCl 
(u'  =  31).  As  for  HgBr,  the  relevant  bands  lie  in  very  con¬ 
gested  regions  of  the  spectra  and  are  likely  blended  to  indis¬ 
tinction  in  spectra  of  "natural”  HgCl  and  Hgl.  The  assign¬ 
ments  can  made  with  confidence  in  our  single  species 


spectra. 

At  present  our  assignments  for  ^****Hg  ”C1  include  39 
bands  spanning  v'  levels  0-9  and  v"  levels  1 1-31.  For  Hgl  we 
have  assigned  50  bands  for  “®Hg  '*’1  and  54  for  ^'“Hg  '”I, 
spanning  v'  =■  0-13  and  v’  =  5-23.  The  assignments  can  be 
extended  provisionally  to  higher  v’  in  both  cases;  however, 
the  bands  in  question  display  anomalous  profiles,  requiring 
care  in  the  estimation  of  band  origins  from  the  measured 
heads.  Consequently,  we  are  continuing  to  work  on  these 
regions.  The  assigned  bands  have  been  least-squares  fitted  to 
the  standard  double  polynomials  in  p{v'  -f  1/2)  and 
piv’  +  1/2)  [see  Eqs.  (1)  and  (2)  in  Ref.  1].  For  Hgl  the  vari¬ 
ance  increases  by  a  factor  of  2  when  the  v’  numbering  is 
altered  by  ±  1  from  our  new  numbering  (which  is  reduced 
by  one  unit  from  Wieland’s),  so  the  new  numbering  is  pre- 


TABLE  I.  Spectroscopic  parameters  (cm  ')  for  the  B-^X  transitions  in 
Hga  and  Hgl.* 


''®Hg”a'' 

“*Hg  '2’I' 

jr. 

23451.6 

24066.4 

<iK) 

191.941 

1 10.850 

-0.4754 

-0.1716 

<i«) 

298.973 

123.053 

-2.1513 

-0.7130 

<> 

1.0112x10-2 

-3.1161x10-2 

c 

-5.0*62X10-* 

2.7567x10  * 

a 

0.37 

0.19 

3>: 

39900* 

38300* 

0: 

8350 

2900 

/t;(A) 

3.02 

3.30 

y?r(A) 

2.42 

2.81 

0.06210 

0.01994 

-2.1X10* 

-4.2x10  ’ 

cr.(«r) 

0.09674 

0.02747 

-  0 

-7.239x10  * 

-  1.954x10  * 

9.006X  10  * 

-6.466x10-’ 

-4.889x10  ’ 

-  1.571x10’ 

*  All  rotational  constants  are  intended  as  guidelines  only,  as  they  are  based 
upon  assumptions  about  the  potential  curves;  see  text. 

'’Vibrational  constants  valid  for  v’  ^  0-9,  v*  =  1 1-31. 

*  Vibrational  constants  valid  for  v’  =  0-13,  u’  =  S-23. 
AssumingdisaociationtoHg* I'S)  +  .Z  “|'S).ThelowestHg*  +  asymp¬ 
tote  lies  17  400cm''  lower  in  HgCl  and  21  800  cm  '  lower  in  Hgl. 


789 


Appl.  Phyt.  Lett.  41(9),  1  November  1982  0003-6951/82/090789-03$01.00 


1 982  Americen  Institute  of  Physics 


789 


TABLE  II.  Isotopicp  bcton  and  abundances  for  isotopic  species  in  HgO  and  Hgl.* 


Hca 

Abundance 

P 

Hgl 

Abundance 

P 

198, 3S 

7.6% 

1.000  732 

198, 127 

10.0% 

1.001  %1 

199, 3S 

12.7 

1.000  374 

199, 127 

16.8 

1.000973 

200,33 

17.4 

1.000000 

20ai27 

23.1 

1.000000 

201,33 

10.0 

0.999  629 

201, 127 

13.2 

0.999032 

202,33 

22.3 

0.999  262 

202, 127 

29.8 

0.998074 

204,33 

3.2 

0.998  338 

204, 127 

6.9 

0.996  181 

198,37 

2.4 

0.977  308 

199,37 

4.1 

0.977  121 

20a37 

3.7 

0.976738 

201,37 

3.2 

0.976  338 

202,37 

7.3 

0.973  982 

204,37 

1.7 

0.973  241 

•  Reference  rntdecubs  are  “®Hg  ”a  and  “Tig  '”1. 


ferred.  For  HgCl  we  have  fitted  only  the  bands  measured  for 
*®“Hg  ”a.  However,  we  have  verified  Wieland’s  numbering 
by  measuring  several  prominent  bands  in  the  v'  =  0  progres¬ 
sion  for  ^‘’’’Hg  ^^Cl.  In  this  case  the  effect  ofthe  change  in  the 
Cl  mass  is  substantial,  so  that  a  change  in  the  v“  numbering 
represents  a  ~  3  cm~ '  shift  in  the  band  positions.  The  agree¬ 
ment  is  within  1  cm~ ',  in  confirmation  of  Wieland’s  num¬ 
bering. 

Minimum  variance  was  achieved  in  the  least-squares 
fits  using  two  upper  state  parameters  and  four  lower  state 
parameters  for  both  molecules.  The  results  are  presented  in 
Table  I.  The  reader  is  cautioned  not  to  trust  these  constants 
outside  the  v'  and  v'  regions  spanned  by  the  assignments. 
For  example,  the  minimum  of  the  .y  state  remains  uncertain 
by  about  20  cm~ '  in  HgG,  relative  to  the  lowest  assigned  v* 
level.  However,  within  the  sampled  regions  these  constants 
should  permit  calculation  of  band  positions  reliable  to  about 
the  standard  deviations  in  the  fits  (0.37  cm~ '  for  HgO,  0.19 
cm~'  for  Hgl).  This  reliability  should  extend  to  the  other 
isotopic  molecules,  for  which  band  positions  can  be  calculat¬ 
ed  by  substituting  the  appropriate  isotopic  p  values  in  the 
polynomials.  These/?  values  are  summarized  in  Table  II.  For 
reference  we  note  that  the  isotopic  shifts  in  the  region  (ff 
strong  emission  are  about  1.1  cm~ '  per  unit  Hg  mass  change 
for  Hgl,  and  about  1.3  cm~ '  for  HgCl. 

We  have  attempted  to  estimate  the  dissociation  energies 
for  the  JIT  states  of  b^  molecules  using  long-range  theory  in 
the  manner  employed  by  Wilcomb  and  Bernstein.^  How¬ 
ever,  in  both  cases  we  have  encountered  the  same  problem 
met  in  the  work  on  HgBr:  for  our  highest  assigned  levels  the 
absolute  sh)pes  in  the  appropriate  long-range  plots  are  al¬ 
ready  greater  than  the  theoretical  limiting  slopes.  Conse¬ 
quently,  we  can  presently  give  only  rough  upper  bounds  on 
.^oO^28S0cm~ '  and  8200 cm~'  for  Hgl  and  HgCl,  respec¬ 
tively.  These  values  are  300-400  cm~ '  lower  than  estimated 
by  Wilcomb  and  Bernstein  and  in  fact  are  close  to  the  origi- 
md  estimates  of  Wieland. 

Franck-Condon  calculations  corroborate  our  vibra¬ 
tional  assignments  and  indicate  that  R,'  is  larger  than  R/ 
by  0.60  k  in  HgCl  and  0.49  A  in  Hgl.  In  these  calculations 
we  Uxdt  the  approach  used  in  the  work  on  HgBr:  we  fixed  the 
X  curves  at  the  values  used  by  Cheung  and  Cool*  and 
varied  the  intemuclear  distaace  in  the  B  states.  The  JT  curves 


were  approximated  as  Morse/RKR  curves,^  and  the  B 
curves  as  Morse  curves.  The  resulting  R, '  values  are  thought 
to  be  reliable  within  0.01  A,  relative  to  the  X  curves  in  the 
region  for  strong  emission.  For  Hgl  our  R/  value  is  0.03  A 
larger  than  Cheung  and  Cool’s.  However,  this  distance  is 
entirely  due  to  a  shift  in  the  attractive  branch  of  the  X  curve 
in  the  Franck-Condon  region  of  strong  emission;  and  with 
the  change  in  the  v’  numbering,  our  Franck-Condon  distri¬ 
butions  are  in  agreement  with  theirs.  For  HgCl  our  '  value 
is  0.08  A  larger  than  Cheung  and  Cool’s.  Most  of  this  differ¬ 
ence  is  significant,  and  it  results  in  a  shift  of  the  Franck- 
Condon  distributions  upward  by  two  v"  units  (e.g.,  the  FC 
gap  for  v'  =  1  occurs  near  v"  =  19  instead  of  v'  =  17,  as 
indicated  in  Table  8  of  Ref.  8). 

The  R,  values  are  known  in  only  a  relative  sense,  as  we 
have  not  yet  analyzed  rotational  structure.  However,  pre¬ 
liminary  examination  of  several  strong  bands  in  the  v'  =  0 
progression  for  HgCl  indicates  that  R/  is  ~3.00  A  in  this 
molecule.  This  value  is  slightly  smaller  than  the  value  in 
Table  I  but  larger  than  Wadt’s'*  theoretical  estimate  (2.93 
A).  For  H  Ar  the  experimental  estimate' '  is  3.06  A,  which  is 
only  0.02  A  larger  than  the  theoretical  value.  Although  the 
rotational  constants  for  Hgl  and  HgCl  have  not  yet  been 
determined  experimentally,  we  include  as  guidelines  in  Ta¬ 
ble  I  the  parameters  calculated  for  the  AT  curves  employed  in 
the  Franck-Condon  calculations. 

Lasing  has  been  reported  for  Hgl  at  a  number  of  wave¬ 
lengths  between  4414  and  4430  A,'^'*  and  for  HgCl  in  the 
3334-3638-A  region.  The  laser  features  generally  cor¬ 

respond  to  the  most  intense  features  in  the  spontaneous 
emission  spectra  and  probably  involve  multipy  overlapped 
rovibronic  transitions  in  the  several  isotopic  molecules  of 
significance  (see  Table  II).  The  prominent  v'-v’  bands  in  the 
lasing  region  for  Hgl  are  0-1 3, 0-14, 1-13, 1-16,2-17,2-18, 
3-18,  3-19, 4-19, 4-20,  3-21,  and  3-23.  The  strongest  fea¬ 
tures  in  the  HgCl  laser  spectrum  occur  near  3380  A  and 
undoubtedly  involve  appreciable  contributions  from  transi¬ 
tions  in  the  nearly  coincident  0-22  and  1-23  bands,  as  has 
been  noted  previously.  The  2-24  and  3-23  bands  could  also 
contribute  significantly,  as  they  lie  very  close  to  the  other 
two  and  have  large  FCFs. 

In  a  recent  paper,  Kvasnik  and  King'^  have  measured 
and  assigned  41  features  in  the  HgCl  laser  spectrum.  While 
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many  of  the  listed  features  likely  involve  some  of  the  indicat¬ 
ed  bands,  we  think  it  is  unwise  to  make  specific  v'-v*  assign¬ 
ments  for  most  of  these  lines,  as  they  probably  involve  cir¬ 
cumstantial  overlap  of  rotational  lines  in  several  bands  of  the 
various  isotopic  molecules.  Furthermore,  Kvasnik  and 
King’s  high-v*  assignments  are  inconsistent  with  our  reana¬ 
lysis  of  this  transition.  These  authors  also  concluded  that  the 
HgCl  emission  in  the  SS40-S730-A  region  includes  a  signifi¬ 
cant  broadband  contribution,  attributed  to  a  bound-free 
transition.  We  think  it  unlikely  that  such  a  bound-free  tran¬ 
sition  could  involve  any  of  the  known  states  of  HgCl  for  the 
following  reasons.  ( 1 )  The  A  transition  should  lie  about 

4000  cm-'  to  the  red  of  J?  — JT.  (2)  The  C-^A  uidD^A 
systems  may  occur  in  this  region;  however,  if  these  systems 
are  present,  the  C— »  JT  and  D-*X  systems  should  occur 
strongly  in  the  UV.  To  our  knowledge  strong  emission  in  the 
latter  systems  has  not  been  reported  for  typical  laser  excita¬ 
tion  conditions.  (3)  The  continuum  of  the  B—^X  transition 
falls  precisely  in  this  region;  however,  the  Franck-Condon 
properties  of  this  system  prohibit  significant  continuous 
emission  from  v’  levels  less  than  1 1.  At  a  typical  operating 
temperature  of  ISO  *C  the  Boltzmann  factor  for  the  sum  of 
all  levels  greater  than  10  is  about  10-^.  Thus  if  vibrational 
thermalization  is  appreciable,  tbeB—t^X  continuum  must  be 
negligible.  While  we  cannot  rule  out  broadband  emission 


from  other  species  in  Kvasnik  and  King’s  laser  system,  we 
think  it  possible  that  their  continuum  is  simply  the  quasicon¬ 
tinuum  of  densely  overlapped  lines  in  the  B—^X  discrete 
emission.  Even  in  our  high-resolution,  single-isotope  spectra 
the  emission  appears  almost  continuous  in  some  regions. 

This  work  was  supported  by  the  Office  of  Naval  Re¬ 
search. 
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DIract  fitting  of  spectroscopic  data  to  near-dissociation 
expansions:  Brj(0' -»>»'),  and  XeCi(B->X  and 

O^X) 

Joel  Tellinghuisen 

Department  of  Chemistry.  Vanderbilt  University,  Nashville.  Tennessee  37235 
{Received  13  October  1982;  accepted  5  November  1982) 

The  utility  of  near-dissociation  expansions  (NDEs)  for  diatomic  vibrational  energies  is  tested  through  least- 
squares  fitting  of  vibrational  bandhead  data  for  selected  electronic  transitions  in  Ij,  Brj,  and  XeCl.  In  these 
test  cases,  the  NDEs  show  efficiency  comparable  to  that  for  the  traditional  polynomials  in  {v  +  1/2).  For  data 
sets  which  span  an  intermediate  range  of  v  levels,  the  NDEs  are  clearly  superior  to  the  polynomials  for 
extrapolating  to  higher  v,  and  of  comparable  reliability  in  extrapolation  to  low  v.  However,  in  their  approach 
to  dissociation,  the  NDEs  can  show  unphysical  behavior,  the  correction  of  which  requires  further  constraints 
on  the  form  of  the  NDE.  Also,  for  lengthy  extrapolations  the  NDEs  may  yield  optimistically  precise  yet 
erroneous  values  for  the  dissociation  energy.  The  present  best  estimates  of  the  dissociation  energies  (£»,)  for 
the  lower  sutes  involved  in  these  calculations  are  2506.0(3)  cm  '  for  Ij(/f  '),  2828(8)  cm  '  for  Brj(/f  ')'  Md 
281  l(7|cm  forXeCl(Jf). 


INTRODUCTION 

Prior  to  1970  the  standard  procedure  for  estimating 
iatomic  dissociation  energies  from  spectroscopic  data 
ras  the  Birge-Siponer  (BS)  extrapolation,*'*  in  which  a 
lot  of  vs  f  is  extended  linearly  to  intercept  the 

axis.  The  area  under  the  extrapolated  curve  is  then 
dded  to  the  energy  of  the  highest  observed  level  to 
leld  an  estimate  of  D,.  A  linear  BS  extrapolation  is 
igorously  correct  for  one  of  the  simplest  “realistic” 
iatomic  potentials— the  Morse  curve.  However,  “real” 
iatomics  do  not  follow  Morse  behavior  near  dissocia- 
lon;  consequently,  the  BS  plot  is  inherently  curved, 
ince  the  BS  method  contains  no  prescription  for  esti- 
lating  the  curvature,  the  resulting  estimates  ofli«  may 
e  very  uncertain. 

To  remedy  this  deficiency,  LeRoy  and  Bernstein*’^ 
sed  semiclassical  theory  to  develop  new  extrapolation 
lethods  based  on  the  asymptotically  limiting  inverse 
ower  potential , 

V(R)»U-C,/R"  .  (1) 

(uch  of  the  early  work  with  these  methods  focused  on 
fie  determination  of  D  and  C,  for  states  for  which  spec- 
roscopic  data  were  available  close  to  the  dissociation 
imit.*  In  such  applications  long-range  methods  yielded 
ery  precise  estimates  of  J),  and  C„  values  which  were 
iostly  in  fair  agreement  with  theoretical  values  calcu- 
Ued  for  the  interacting  atoms. 

As  it  turns  out,  the  long-range  limiting  behavior  is 
ot  very  sensitive  to  errors  in  C„.  Taking  advantage 
f  this  fact,  Tellinghuisen  el  al.*  and  Wilcomb  and 
lernstein*  employed  a  variant  method,  in  which  C„  was 
ixed  at  its  theoretical  value,  to  estimate  Xi’s  from  data 
^hich  required  long  extrapolations.  In  a  recent  paper 
eRoy  and  Lam*  have  taken  this  idea  further,  using  non- 
inear  least-squares  fits  to  expressions  having  the  cor- 
ect  limiting  behavior  at  dissociation,  with  values  of 
and  C,  again  fixed  by  theory. 

In  the  present  work  I  have  expanded  on  the  method  of 


LeRoy  and  Lam®  in  a  test  of  direct  fitting  of  band- 
head  data  to  near-dissociation  expansions.  Since  all 
spectroscopic  data  represent  energy  differences,  direct 
fits  of  raw  data  to  expressions  which  include  the  energy 
levels  for  all  involved  states  are  preferred  over  methods 
which  attempt  to  isolate  the  dependences  on  the  different 
states  through  preliminary  manipulation  of  the  raw  data. 
The  latter,  which  include  combination  difference  and 
term  value  methods.  Introduce  bias  and  correlation 
error  Into  the  final  results.®’*®  In  the  treatment  of 
BeAr*  by  LeRoy  and  Lam,®  the  fitted  energy  levels  were 
already  a  step  removed  from  the  raw  data  of  Subbaram 
et  al.  **  In  the  present  case  I  have  fitted  bandhead  data 
of  electronic  emission  transitions  in  Ij,  Br^,  and  XeCl. 
These  transitions  all  involve  low  v'  levels,  so  the  tra¬ 
ditional  polynomial  in  (v'  +i)  is  used  to  represent  G'„  in 
each  case.  For  XeCl  the  data  sample  a  large  fraction 
of  the  ground  state;  for  I,  and  Br^  an  intermediate  range 
of  tf'  is  sampled.  In  the  calculations  I  have  focused 
my  attention  on  the  following  question:  How  do  (1)  the 
efficiency  and  (2)  the  extrapolating  ability  of  the  near- 
dissociation  expansions  comjpare  with  those  of  the  con¬ 
ventional  polynomials?  The  results  are  generally  fav¬ 
orable  for  the  long-range  method  on  both  counts;  how¬ 
ever,  some  limitations  are  noted. 

II.  THEORY 

Near  the  dissociation  limit  the  vibrational  energy  fol¬ 
lows  the  expression,  ® 

,  (2) 

where  Vp  is  the  effective  vibrational  quantum  number  at 
dissociation,  n  is  the  theoretically  appropriate  inverse 
power  of  6),  and  the  constant  X„  contains  the  de¬ 
pendence  on  C„  and  the  reduced  mass  4, 

.  (3) 

Vaiues  for  the  constant  X„  for  various  n  may  be  found 
in  Ref.  8. 

A  different  version  of  Eq.  (2)  is 
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dGjdv  =  A,(5),  -  GJ ,  (4) 

where*  Jf,  =[2n/(w  -  A  modified  form  of 

Eq.  (4), 

idGjdv)'^'^^^  =/:?''*’**>(»,  -  G.) .  (5) 

was  employed  in  the  graphical  extrapolation  methods  of 
Refs.  6  and  7. 

LeRoy  and  Lam*  suggested  the  fitting  of  G,  values  to 
an  empirical  extension  of  Eq.  (2), 

G,=a),  ,  (6) 

in  which  the  correct  limiting  behavior  is  assured  by 
designing  the  function  Fiv)  to  go  to  unity  as  v^Vd‘  The 
two  functions  examined  in  Ref.  8  were 

=  l  +  ^  -«)'  (7) 

and 

["•  'Ih/h-S)  /o^ 

1  +  ^6|(«d -wl'J 

In  Eqs.  (7)  and  (8)  the  {aj  and  {&,}  are  purely  empirical 
quantities  to  be  determined  from  a  fit  of  the  data.  In 
further  work  on  rotational  and  centrifugal  distortion 
constants,  Tromp  and  LeRoy‘*  have  employed  exponen¬ 
tials  for  F{v).  In  the  present  work  1  have  used  only  the 
J^s  of  Eqs.  (7)  and  (8)  or  simple  modifications  thereof. 

For  the  four  band  systems  investigated  here,  previous 
work  has  shown  that  two  vibrational  constants  are  ade¬ 
quate  to  represent  the  v'  levels.  Thus  the  data  are 
fitted  to  expressions  of  the  form , 

v{v\v")^T^G[*X„(vd-v")^'^-*^F(v'>)  ,  (9) 

where  G' =  +5)  -  «, +i)*.  Note  that  the  param¬ 

eter  T  represents  the  energy  difference  between  the 
minimum  of  the  upper  state  and  the  dissociation  limit 
of  the  lower  state;  in  conventional  polynomial  fitting 
the  corresponding  term  is  . 

The  nonlinear  fits  were  carried  out  using  standard 
methods.*®'”  In  particular  the  bookkeeping  was  great¬ 
ly  facilitated  by  use  of  the  coefficient  matrix  U,  which 
contains  m  columns  (one  for  each  adjustable  parame¬ 
ter)  and  n  rows  (one  for  each  experimental  point) .  ‘®  In 
these  nonlinear  fits  the  elements  of  the  fth  row  of  this 
matrix  are  the  partial  derivatives  of  the  fit  function  with 
respect  to  the  parameters,  evaluated  at  the  ith  point  us¬ 
ing  the  current  estimates  of  the  parameters.  The  cor¬ 
rections  to  the  parameters  are  then  calculated  straight¬ 
forwardly  upon  inverting  the  matrix  U*'  U;  and  the  pro¬ 
cedure  is  repeated  until  convergence  is  obtained.  ** 
Variable  weights  are  readily  accommodated. 

Also  of  interest  are  the  errors  in  the  function  and  in 
other,  derived  functions  of  the  parameters.  The  errors 
are  calculated  from  the  variance-covariance  matrix 
[which  is  proportional  to  (U®’  U)'*],  using  methods  such 
as  those  described  in  Refs.  13  and  14.  The  quantities 
of  particular  interest  are 

Ar,  =  i/(-0.5, -0.5)  (10) 


and 

t>y  =  v{v' , -O.i) -v(v' ,Vd)  .  (11) 

III.  RESULTS  AND  DISCUSSION 
A.  Ij(D'-A') 

I  have  recently  reported  a  detailed  vibrational  analysis 
of  the  D'(2g)  -A'(2tt)  system  (~  3400  A)  of  Ig,  **  with 
~  250  assigned  bands  for  **’lg  and  ***Ij,  spanning  c'  lev¬ 
els  0-15  and  r"  levels  4-30.  By  identifying  the  D' 
state  as  the  a  state  of  King  et  al. ,  *®  I  was  able  to  esti¬ 
mate  that  S>,  for  the  A!  state  is  very  close  to  2500  cm'*. 
In  more  recent  work  Koffend  et  al.^''  have  achieved  a 
complete  vibrational  and  rotational  analysis  of  the  A' 
state.  Their  work  corroborates  my  assignment  of  the 
emission  spectrum  and  yields  precise  information  for  v" 
levels  0-59.  Thus  this  system  makes  an  interesting 
test  for  fitting  to  near -dissociation  expansions  (NDEs). 

Data  for  multiple  isotopic  molecules  can  be  readily  ac¬ 
commodated  in  fits  to  Eq.  (9),  just  as  they  can  in  con¬ 
ventional  polynomial  fitting.  However,  the  vibrational 
numbering  of  the  A'  state  is  now  certain,  so  data  for 
multiple  isotopes  are  not  needed.  Hence  I  have  fitted 
only  the  **%  bands  in  Table  I  of  Ref,  15,  with  their  as¬ 
sociated  weights. 

For  the  A'  state  «  =  5.  For  C5  I  have  used  the  value 
Cj  =2.64x10*  cm'*  A*,  given  by  Mulllken**  and  used  in 
Ref,  1 5.  To  obtain  initial  estimates  of  the  parameters 
I  used  linearized  versions  of  Eqs.  (7)-(9).  In  this  re¬ 
gard  it  is  worth  noting  that  F,,„(t')  and  variants  thereof, 
with  selected  powers  of  (t^^  -  v)  omitted,  are  particularly 
convenient,  since  with  Vp  fixed,  the  fits  to  Eq.  (9)  be¬ 
come  linear  (which  means  that  convergence  occurs  in 
one  pass).  For  comparison  with  the  NDE  fits,  I  have 
also  redone  the  polynomial  fits  for  the  *®*l2  data  alone. 

The  results  of  these  calculations  are  summarized  in 
Table  I.  In  the  polynomial  fits  three  constants  were 
clearly  inadequate  for  the  A'  state,  and  the  most  effi¬ 
cient  representation  was  the  P*,,  fit  (i.e.,  two  upper- 
and  four  lower-state  parameters,  plus  AT,).  In  the 
NDE  fits  the  number  of  parameters  ranged  from  six 
to  eight;  in  each  case  two  vibrational  constants  were 
used  for  the  upper  state.  The  six-parameter  NDE  fits 
actually  outperformed  the  six-parameter  polynomial 
fit.  However,  the  variances  for  the  seven-parameter 
NDE  fits  were  noticeably  higher  than  that  for  the 
fit.  In  fact  I  was  unable  to  make  the  nonlinear  F,^^ 
and  fits  converge.  However,  by  treating  Vg  as  a 
known,  then  varying  it  externally,  I  was  able  to  find  a 
minimum  variance  with  respect  to  Vp.  To  calculate  the 
errors  in  the  nonlinear  fits,  I  then  simply  froze  the  pa¬ 
rameters  at  the  values  determined  from  optimization 
with  respect  to  Vp.  The  instability  in  the  seven-param¬ 
eter  NDE  fits  was  manifested  as  large  errors  in  and  cor¬ 
relation  between  T  and  Vp.  The  F„,4  fit  also  failed  to 
converge,  and  in  this  case  the  minimum  variance  with 
respect  to  a  "known”  Vp  occurred  at  the  unrealistically 
large  value  228,  at  which  point  T  went  to  zero. 
Somewhat  surprisingly  the  variance  in  this  fit  matched 
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TABLE  I.  Summary  of  various  least-squares  fits  for  the  I2  D'—A'  system.* 


Fit” 

(tr»)“ 

AT," 

< 

■n"‘* 

"c 

^2.3 

0.669 

30  367.2(1.7) 

103.99(9) 

0.210(8) 

... 

... 

P2,4 

0.102 

340.1(1.2) 

4. 00(4) 

0.210(3) 

... 

... 

P2.5 

0.101 

337.4(2.7) 

3.99(4) 

0.209(3) 

... 

... 

P3.4 

0.102 

340.0(1.3) 

4.06(7) 

0.224(14) 

... 

... 

^3,5 

0.102 

337.8(2.7) 

4.04(7) 

0.222(14) 

... 

... 

^’..2 

0.462 

357.7(1.9) 

4.05(8) 

0.213(6) 

2373(6) 

63.4(2) 

^Oi3 

0.262 

348.0(5.1) 

4.02(6) 

0.211(5) 

2516(890) 

82(142) 

Fa.i 

0.101 

335.5(1.4) 

3.98(4) 

0. 209(3) 

... 

... 

0.472 

358.4(1.9) 

4.05(8) 

0.213(6) 

2372(6) 

63.3(3) 

^b,3 

0.245 

346.2(5.6) 

4.02(6) 

0.211(5) 

2538(2598) 

90(905) 

Fa.i 

0.483 

358.3(2.0) 

4.05(8) 

0.214(6) 

2364(6) 

61.6(2) 

Fa.3 

0.229 

346.1(4.7) 

4.02(6) 

0.211(5) 

2545(793) 

81(97) 

*AU  quantities  in  cm"'  except  a’icm"®)  and  vj,  (dimensionless). 

refers  to  polynomial  fit  having  i  upper-state  parameters,  J  lower-state  parameters,  plus 

AT,. 

“Defined  as  [(rM'j6i’j)/(Ztrj)l  x  [n/(n  —  m)l,  where  nq  and  fivi  are  the  weight  and  residual  of  the 
ith  point,  «  is  the  number  of  data  points,  and  w  is  the  number  of  parameters. 

“’current  best  estimates,  A7i  =  30 347.  0cm"',  l)',' --- 2506. 0  cm"*;  see  Ref.  18. 

•Fits  for  altered  Cj  value;  see  the  text. 


that  of  the  best  polynomial  fitl  However,  the  NDE  pa¬ 
rameters  were  physically  meaningless. 

In  all  fits  the  vibrational  parameters  for  the  D'  state 
were  well  determined  and  mutually  consistent.  Further¬ 
more  the  estimates  of  AT,  from  the  NDE  fits  were 
close  to  those  from  the  polynomials  and  showed  com¬ 
parable  scatter  and  precision.  In  fact  the  values  from 
the  seven-parameter  fits  are  very  close  to  my  current 
best  estimate,  30347. 0  cm"’,  obtained  by  correcting  my 
results*®  in  accord  with  the  better  definition  of  the  A' 
state  in  Ref.  17.*® 

Some  of  these  results  are  displayed  in  Fig.  1,  which 
is  the  plot  suggested  by  Eq.  (15)  and  used  previously  in 
my  laboratory  for  graphical  estimation  of  31,.  I  have 
included  also  the  points  calculated  from  the  tabulated 
energy  levels  given  in  Ref.  17,  using  dCv/dv^{G^i 
-  G„.,)/2 .  The  latter  points  diverge  from  ray  polynomial 
results  above  r"  =25,  suggesting  that  my  assignments 
for  v"  >25  may  be  in  error.  When  I  dropped  the  ques¬ 
tionable  bands  from  the  data  set,  the  variances  in  the 
nonlinear  fit  to  ^.,3  indeed  dropped  by  a  factor  of  2; 
however  the  convergence  problem  persisted.  From  the 
data  of  Ref.  17  it  is  clear  that  the  jog  in  the  curve  near 
v  =  23  is  real,  even  though  my  data  may  overemphasize 
it.  This  points  up  a  limitation  inherent  in  any  attempt 
to  estimate  31,  from  data  which  are  far  from  dissociation. 
There  is  clearly  no  way  for  any  least-squares  tit  to  cor¬ 
rectly  predict  this  type  of  behavior  from  data  which  do 
not  sample  the  relevant  region.  In  fact  below  t;  =  22  the 
data  in  Fig.  1  appear  to  follow  the  calculated  limiting 
slope  very  nicely;  lacking  data  above  v  =  22,  one  would 
estimate  a  1>,  value  which  is  much  below  the  true  value 
and  optimistically  precise. 

Another  intersting  point  about  the  NDE  results  in  Fig. 

1  is  that  even  though  they  are  constrained  to  have  the 
correct  limiting  behavior,  the  seven-parameter  fits 
appear  not  to  approach  dissociation  with  the  correct 
slope.  In  fact  they  do,  but  only  very  near  wj,.  This 


point  is  discussed  further  below,  in  conjunction  with 
the  Brj  calculations. 

To  eliminate  model  dependence  in  the  estimates  of 
X),  and  Vj,,  LeRoy  and  Lam®  took  weighted  averages  of 
the  various  estimates,  with  weights  proportional  to  the 
reciprocal  variances  of  the  individual  estimates.  In  the 
present  case  this  procedure  would  mean  that  the  values 
are  almost  entirely  determined  by  the  six-parameter 


V 


FIG.  1.  Plots  oig  (units  cm"'®'’)  vs  G„  for  the  A’(2u)  stale  of 
I2.  The  large  points  are  results  from  the  ^2,1  polynomial  fit; 
the  small  points  are  from  results  of  Ref.  17.  Also  shown  are 
the  curves  calculated  from  the  results  of  two  other  polynomial 
fits  and  five  near-dissociation  fits.  The  limiting  theoretical 
behavior  is  indicated  by  the  dashed  line  (T) .  All  calculated  g 
curves  have  their  G,  scales  adjusted  to  be  coincident  at  20. 
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FIG.  2.  Plots  of  Y  =  AT,-G„  and  Its  error  oy  as  (unctions  of 
V  for  the  A’  state  of  Ij.  The  solid  curves  represent  results  of 
the  ^2,4  polynomial  fit,  dashed  curves  the  near-dissocia¬ 
tion  fit.  Note  the  logarithmic  scale  for  ay. 


fits,  since  and  Vg  are  so  ill-determined  in  the  seven- 
and  eight-parameter  fits.  The  result  would  be  an  ap¬ 
parently  very  precise  value  of  2372. 6  cm*‘  for  which 
is  130  cm'*  (or  about  i  the  extrapolation  energy)  below 
the  correct  value.  It  is  clearly  necessary  to  more  fully 
sample  the  model  dependence.  A  way  of  doing  so  with 
the  framework  of  the  F,  and  forms  is  to  drop  various 
powers  of  {vo  -v).  This  method  is  explored  further  in 
the  Br,  work  discussed  below. 

One  troubling  feature  associated  with  model  depen¬ 
dence  is  the  observation  that  bad  fits  having  few  param¬ 
eters  can  often  yield  apparently  more  precise  estimates 
than  ^oo4f  fits  having  more  parameters.  For  example 
the  six-parameter  fits  in  Table  I  all  yield  systematically 
high  values  of  AT,,  with  precision  comparable  to  or 
better  than  that  of  the  seven-,  eight-,  and  nine-param¬ 
eter  fits.  A  second  case  of  "model  deception!’  is  illus¬ 
trated  in  Fig.  2,  which  shows  the  quantity  {AT,  -G”) 
and  its  error  band  for  the  P*,,  and  F„,  fits.  The  error 
for  the  latter  is  larger  everywhere,  even  though  the  estima¬ 
tion  of  levels  between  the  highestobservedv"  level  and  vp 
must  surely  be  better  for  the  NDE.  To  avoid  such 
problems  it  may  be  wise  to  confine  the  model  weighting 
to  fits  having  form  and  quality  comparable  to  that  of  the 
best  fit. 

In  one  other  test  calculation  using  the  I|  data,  I  inves¬ 
tigated  the  effect,  of  changing  the  C5  value.  For  the 
B{0^  *n)  state  of  the  best  experimental  estimate*"  of 
Cj  is  a  factor  of  1. 56  smaller  than  the  value  given  in 
Ref.  19.  Reduction  of  the  C(  for  the  A'  state  by  the 
same  factor  yields  a  value  of  1. 7x  lo*  cm'*  A',  which 


increases  X,  by  34%  and  the  limiting  slope  in  Fig.  1  by 
13%.  The  NDE  fits  conducted  with  this  X,  were  com¬ 
parable  in  quality  to  the  previous  fits,  and  the  conver¬ 
gence  problem  persisted  for  more  than  six  parameters. 
In  keeping  with  expectations  the  estimates  of  changed 
by  only  9  cm'*  for  the  F,,^  fit  and  29  cm'*  for  the  (very 
uncertain)  F,,^  fit. 

B.  Br2(D'->-A') 

A  detailed  analysis  of  the  D'  ^A'  system  in  Brj  has 
been  published  recently.  **  The  final  recommended  con¬ 
stants  for  this  system  were  obtained  from  a  global  fit 
of  all  assigned  rotational  lines  and  band  heads  for  the 
(79,79),  (79,81),  and  (81,  81)  isotopic  Brj  molecules, 
and  spanned  v'  levels  0-6  and  v"  levels  5-21.  The 
highest  observed  v”  levels  appeared  to  follow  the  cal¬ 
culated  limiting  behavior  in  the  plot  equivalent  to  Fig.  1, 
and  a  linear  extrapolation  gave  I>,  =  2835  cm'*  and  vp 
=  53. 9.  To  avoid  the  additional  complications  of  includ¬ 
ing  rotational  structure  and  multiple  isotopes,  I  have 
used  in  the  present  calculations  just  the  red -degraded 
bandheads  for  **Br2  listed  in  Table  1  of  Ref.  21.  These 
bands  sample  the  somewhat  smaller  v"  range,  5-16. 

For  C5 1  have  used  the  same  value  used  in  Ref.  21 
(1.39x10*  cm'*  A*). 

These  calculations  showed  behavior  similar  to  that 
described  for  the  I*  data.  Four  vibrational  parameters 
were  again  needed  for  an  adequate  representation  of  G” 
in  the  polynomial  fitting,  and  two  sufficed  for  the  upper 
state.  Additional  parameters  gave  no  significant  im¬ 
provement.  In  the  NDE  fits,  five  parameters  were 
clearly  inadequate,  and  the  fits  containing  more  than 
six  again  failed  to  converge.  The  six-parameter  fits 
gave  variances  only  slightly  higher  than  the  benchmark 
P*.,  fit. 

Figure  3  illustrates  the  results  of  these  calculations. 
The  aforementioned  problem  of  the  incorrect  slopes  is 
even  more  evident  here.  Figure  4  shows  that  in  the  re¬ 
gion  very  close  to  Vp  the  slopes  approach  the  correct 
values,  as  indeed  they  must.  However,  it  is  clear  that 
the  behavior  is  deviating  from  the  limiting  behavior 
much  too  rapidly  with  increasing  {vp  -v),  because  the 
correction  functions  are  “turning  orf’  too  quickly.  To 
produce  a  more  reasonable  extrapolation,  I  eiqierimented 
with  higher  lead  powers  of  {Vp  -  v)  in  F,.  The  curves 
for  P^2  (2, 4)  [i.  e. ,  terms  of  powers  2  and  4  in  {Vp  -  v)] 
in  Figs.  3  and  4  show  that  dropping  the  linear  term  pro¬ 
duces  considerable  improvement.  The  use  of  still 
higher  powers  gave  even  longer  near-linear  extrapola¬ 
tions,  with  the  estimated  3D,  value  approaching  ~  2820 
cm'*  for  the  highest  powers  investigated  [F,,t(7, 8)].  The 
variances  of  all  of  these  six-parameter  fits  having  no 
linear  correction  term  were  only  marginally  higher 
than  that  for  the  seven-parameter  Pg,,  fit. 

With  the  importance  of  choice  of  lead  power  so  evi¬ 
dent  1  returned  to  the  five-parameter  NDE  fits  [i.  e. , 
those  having  only  one  term  in  ( Vp  -  v)^  in  the  correction 
function].  Remarkably  the  variance  decreased  progres¬ 
sively  as  j  was  increased,  until  for  7  =  7  this  five-pa¬ 
rameter  fit  actually  outperformed  the  seven -parameter 
polynomial  fiti  Further  increase  in  j  led  to  immediate 
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FIG.  3.  g  vs  G„  for  the  A'{2u)  state  of  Brj.  The 

dashed  line  (T)  indicating  the  theoretical  limiting  behavior  is 
arbitrarily  drawn  to  intercept  the  axis  at  3000  cm'*.  Other 
curves  are  as  identified  in  Fig.  1  and  are  adjusted  to  have 
coincident  G„(15)  values. 


deterioration  in  the  fit  quality.  As  j  was  increased 
from  one  to  nine,  the  estimate  of  decreased  mono- 
tonically  from  3091.6(6.8)  to  2751.8(2.4)  cm*^; concomi¬ 
tantly,  AT,  decreased  from  39  747. 6(4. 3)  to  35  594. 6(3. 5) 
cm"*.  (Because  of  the  changing  AT,,  the  extrapolated  ener¬ 
gy  from  the  highest  observed  level  actually  varies  only  half 
as  much  as  does  . )  In  keeping  with  the  previous  comments 


C.  XeCKfl-XandD-X) 

For  these  calculations  I  have  used  the  bandhead  mea¬ 
surements  for  *’*Xe  ’*CI  rejported  in  Ref.  22.  These 
data  span  v  =  0-12  in  the  B  state,  0-9  in  D,  and  0-13  in 
X.  Here,  n  =6,  and  I  have  used  the  C,  value  of  8x  10* 
cm"*  A*  estimated  in  Ref.  6.  The  optimal  polynomial 
fit  of  Ref.  22  contained  ten  parameters— T„  u,,  and 
u>,x,  for  each  upper  state,  and  four  vibrational  param¬ 
eters  for  the  X  r^ate.  Graphical  extrapolation  yielded 
an  estimated  X>,  of  281(10)  cm"*. 

In  the  NDE  fits  both  transitions  were  fitted  simulta¬ 
neously,  as  they  were  in  the  polynomial  fits.  For  the 
X  state,  correction  functions  containing  only  one  term 
in  {vg  -  vY  were  not  adequate.  However  all  two-  and 
three-term  functions  (nine-  and  ten-parameter  fits) 
gave  convergence,  with  variances  mostly  within  10%  of 
that  for  the  polynomial  fit.  Since  the  data  set  included 
v”  =  0,  there  was  very  little  spread  (0.3  cm'*)  in  the 
T,  estimates.  As  before,  there  was  complete  consis¬ 
tency  in  the  various  estimates  of  the  B-  andZi-state  vi¬ 
brational  constants.  In  this  case  the  extrapolation  from 
the  highest  observed  level  of  the  X  state  to  a,  is  short, 
and  the  fits  were  only  weakly  sensitive  to  choice  of 
powers  of  {vg  -  v)  in  the  correction  function.  For  ex¬ 
ample  in  the  nine-parameter  fits  to  ]>,  decreased 
from  285. 4(8)  cm'*  to  277. 4(6)  cm"*  as  the  powers  ot 
{vg  -v)  were  increased  from  (1,2)  to  (3,6).  Weighted 
averages  of  results  from  11  fits  yielded  X>,  =  281. 1(7) 
cm'*,  t)o  =  19. 30(10). 

The  original  analysis  of  the  XeCl  B~X  system*  em¬ 
ployed  sources  containing  Xe  and  Cl  in  natural  isotopic 
abundance.  The  spectra  sampled  the  smaller  v"  range 
0-7  and  yielded  (by  graphical  extrapolation)  a  signifl- 


about  problems  with  model  dependence,  it  is  worth  noting 
that  the  error  estimate  on  the  unreasonably  large  first 
D,  value  is  a  factor  of  2  smaller  than  those  on  the  more 
reasonable  values  from  the  six-parameter  fits.  The 
optimal  five-parameter  fit  gave  the  following  results; . 
AT,  =35645.9(1.8)  cm'*,  D,  =  2820(1.0)  cm'*,  vg 
=  48.444(10),  and  0,=  - 1.7821X10'**.  The  AT,  value 
is  15-55  cm'*  smaller  than  the  estimates  from  the  poly¬ 
nomial  fits,  which  in  fact  do  not  pin  this  quantity 
down  very  well. 


As  for  the  Ig  calculations,  all  NDE  fits  having  vari¬ 
ances  within  a  factor  of  2  of  the  minimum  gave  mutually 
consistent  values  for  the  upper  state  vibrational  param¬ 
eters.  Weighted  averages  from  these  same  fits  yielded 
the  following  results:  AT,  =  35  654(5)  cm'*,  X).,  =  2828(8) 
cm'*,  and  Uo  =48.42(7).  Although  these  averages  were 
obtained  from  results  of  14  fits,  they  were  dominated 
by  the  contributions  from  the  three  included  five-pa¬ 
rameter  fits  F,,i(6),  F,,,(7),  andF,,,(8).  TheD,  and 
Vg  values  are  reasonably  close  to  those  obtained  from 
the  full  data  set  by  graphical  extrapolation.  **  Of  course, 
their  validity  rests  on  the  assumption  that  behavior  such 
as  that  observed  for  Ig(A')  does  not  set  in  above  the  high¬ 
est  observed  level.  In  view  of  that  possibility  for  the 
lengthy  extrapolation  (*-400  cm'*  from  v  =  2\),  the  error 
estimates  are  probably  optimistic. 


FIG.  4.  Near-dlssoclatlon  curves  of  Fig.  3  In  the  region  very 
near  the  dissociation  limit.  The  G„  scales  have  been  adjusted 
to  emphasize  the  comparisons. 
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cantly  smaller  estimate  of  3>,  =  255(10)  cm'‘.  On  repro¬ 
cessing  these  data  with  the  present  direct  fitting  method, 
I  obtained  273(4)  cm'*,  which  is  significantly  closer 
to  the  current  best  estimate.  For  comparison  a  refit 
of  the  1078  data,  truncated  to  the  same  bands  included 
in  the  1975  analysis,  yielded  3),  =  278(3)  cm'*,  which  is 
within  one  standard  error  of  the  value  obtained  from  the 
full  data  set. 

IV.  SUMMARY 

From  these  calculations,  I  draw  the  following  con¬ 
clusions: 

(1)  Near -dissociation  expansions  can  represent  vi¬ 
brational  energies  with  efficiency  comparable  to 
that  of  the  customary  polynomials  in  (v-i-i).  This 
result  bears  out  indications  of  previous  studies^* 

(in  which,  however,  no  explicit  NDE -vs -polynomial 
comparisons  were  given).  For  data  sets  lacking  data 
for  low  V  levels,  the  extrapolations  to  v  show  preci¬ 
sion  and  model-dependent  scatter  comparable  to  those 
for  polynomial  fits.  For  extrapolation  to  large  v,  the 
NDE  fits  must  perform  better,  as  they  are  designed  to 
approach  dissociation  correctly.  Since  these  fits  are 
nonlinear,  convergence  problems  may  occur,  especially 
for  higher  fit  dimensions. 

(2)  The  near-dissociation  fits  may  not  yield  very  good 
long-range  behavior,  even  though  the  NDEs  are  con¬ 
strained  to  have  the  correct  n  and  C,  values.  Their 
performance  can  be  improved  by  excluding  low  orders 

of  (vg  ~v)  from  the  correction  functions.  There  is  some 
theoretical  justification  for  this  step,  **  from  considera¬ 
tion  of  higher-order  contributions  to  the  long-range 
potential.  Although  one  could  incorporate  additional 
constraints  in  the  NDE  to  reflect  these  higher -order 
contributions,  such  a  move  is  hardly  warranted  in  the 
present  examples,  where  the  vibrational  levels  sampled 
by  the  data  are  far  from  the  dissociation  limit.  How¬ 
ever,  it  is  advisable  to  at  least  examine  the  NDEs 
graphically,  to  ensure  that  unacceptable  behavior  like 
that  of  some  of  the  curves  in  Figs.  1  and  3  be  avoided. 

(3)  There  is  no  simple  prescription  for  quantitative 
assessment  of  model  dependence,  as  the  range  of  pos¬ 
sible  correction  functions  is  essentially  unlimited.  To 
avoid  some  of  the  problems  of  model  error,  it  may  be 
wise  to  consider  in  the  final  evaluation  only  fits  having 
quality  comparable  to  that  of  the  “best”  fit.  In  the 
present  calculations,  only  the  correction  functions  of 
Eqs.  (7)  and  (8)  have  been  examined.  Rational  poly¬ 
nomials**"*  or  exponentials**  may  prove  better  in  many 
applications. 

(4)  In  the  cases  investigated  here,  the  vibrational  con¬ 
stants  for  the  upper  states  were  consistent  with  the  poly¬ 
nomial  results,  for  all  NDE  fits  having  variances  com¬ 
parable  to  that  of  the  best  fit.  This  result  suggests  that 
one  can,  with  negligible  bias  and  correlation  error,  fit 
to  a  modified  form  of  Eq.  (9),  in  which  G'  is  taken  as 
known.  Howevt  one  should  nof  fit  the  polynomial - 
derived  G"  expression  directly,  as  this  quantity  may 
already  contain  considerable  model  dependence.  This 
point  is  particularly  relevant  for  limited  data  sets,  such 


as  those  for  the  D' -A'  systems  of  Ij  and  Br*,  where  AT, 
and  G"  are  strongly  dependent  on  polynomial  order. 

(5)  The  achievement  of  a  good  NDE  fit  for  a  limited 
data  set  should  never  be  accepted  as  an  excuse  for  not 
seeking  additional  data.  The  reason  is  that  no  fit  can 
be  expected  to  correctly  account  for  behavior  such  as 
that  of  IfiA*)  in  Fig.  1  unless  the  data  sample  the  rele¬ 
vant  regions.  This  type  of  behavior  may  be  more  the 
rule  than  the  exception,  since  it  occurs  far  from  the 
asymptotic  region  where  long-range  theory  should  apply. 
In  fact,  work  in  my  laboratory  has  shown  similar  effects 
for  the  A  state  of  Ij*’*’  and  the  X  states  of  the  mercury 
halides*’**:  In  the  former  case  the  plot  of  Eq.  (5)  dis¬ 
plays  a  slope  which  is  about  half  the  theoretical  slope, 
until  very  close  to  the  limit;  in  the  latter  the  slopes 
exceed  the  limiting  slopes  for  intermediate  v,  then  ap¬ 
proach  the  limit  from  above.  In  each  case  NDE  fits 
to  limited  data  sets  would  yield  systematically  erroneous 
3)«  values.  In  this  light  it  may  be  provident  to  report 
errors  in  these  determinations  with  an  ounce  of  pes¬ 
simism. 
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The  B-X,  C-X,  and  D-X  cranslclons  of  HgCl,  HgBt,  and  Hgl  are 
reanalyzed  using  Tesla  discharge  sources  containing  single  Isotopic 
species  of  these  molecules.  Direct,  simultaneous  least-squares  fits 
of  all  transitions  yield  optimal  vibrational  constants  for  all  four 
states.  Low-resolution  studies  of  the  emission  as  a  function  of 
buffer  gas  pressure  show  effects  of  vibrational  relaxation  In  the  B 
state  and  colllslonal  quenching  of  the  C  and  0  states.  The  broad  B-A 
bands  In  the  red  and  Infrared  are  weaR  but  clearly  present.  In 
support  of  a  recent  report  of  this  transition  In  HgBr.  In  addition 
there  Is  evidence  of  ocher  cranslclons,  previously  unreported  for 
these  molecules.  Including  a  very  weak,  system  near  2200  A  In  Hgl, 
which  shows  fine  red-degraded  vibrational  band  structure, 

INTRODUCTION 


Although  Che  mercury  halide  lasers  are  of  considerable  current 
Interest,  surprisingly  little  spectroscopic  work  has  been  done  on  the 
2  +  2  + 

lasing  B(  £  )  -  X<  z  )  cranslclons  in  Che  HgX  molecules.  The  main 
source  of  Infonsaclon  on  these  systas  has  been  a  series  of  papers 

published  by  Wleland  over  20  years  ago,^”^  Most  of  Wleland's  work 
Involved  sources  containing  Che  mercury  halides  In  natural  Isotopic 
abundance.  To  refine  Che  spectroscopic  characterization  of  the  HgX 
molecules,  my  group  has  been  reanalyzing  their  emission  spectra, 
using  sources  prepared  from  single  Isotopes  of  Hg  and  the  halogens. 
Preliminary  reports  of  our  work  on  the  B-X  systems  of  HgCl,  HgBr,  and 

Hgl  have  been  published.  ’  In  this  paper  I  discuss  our  continuing 


2  2 

work  on  these  systems  and  new  work  on  the  C(  II,,,)-X,  D(  Hi/o)"*,  and 

5  A/ j/ 4 

B-A(  u)  systems.  In  addition  we  see  evidence  of  weaker  transitions, 
not  previously  reported. 


HlGH-RESOLbTION  STUDIES 


Emission  spectra  have  been  photographed  at  high  resolution  from 
sources  containing  the  Individual  Isotopic  species,  ^^^Hg^^Cl, 


200Hg37ci,  200Hg'^r,  ^OO^^S^Br,  200„gl27^^ 


excitation  by  means  of  a  Tesla  discharge.^  Details  of  Che 
preparation  and  operation  of  the  sources  have  been  given  In  the 

4  3 

preliminary  reports.  ’  The  sources  show  optimal  emission  in  the 
visible  B-X  systems  for  HgX^  pressures  around  1  torr  and  Ar  buffer 
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The  8(^2^)  X(22*)  Transition  (4050-4500  A)  in  Hgl 

K.  S.  ViSWANATHAN  AND  Joel  Tellinohuisen 
Department  of  Chemistry,  yanderbilt  University,  Nashville.  Tennessee  J72JS 

The  B  X  band  system  (4030-4 $00  A)  of  Hgl  is  photographed  and  vibrationally  analyzed 
for  the  isotopically  pure  species  ’™Hg‘'’l  and  *^fg'”l.  The  assigned  bands  span  tf  levels  5-26 
and  1/  levels  0- 1 3.  The  least-squares  analysis  indicates  that  the  previously  accepted  numbering 
for  this  system  is  one  unit  too  high.  Band-profile  simulations  and  Eranck-Condon  calculations 
indicate  that  the  intemuclear  separation  i$  0-49  A.  The  ground-state  dissociation  energy 

(2),)  is  estimated  to  be  2750  ±  80  cm~'.  Spin  splitting  is  found  to  contribute  significantly  to  the 
band  structure. 

INTRODUCTION 

The  B  —*X  emission  systems  of  the  diatomic  mercury  halides  (HgX)  have  become 
a  subject  of  renewed  interest  in  recent  years,  following  the  discovery  of  strong  lasing 
on  these  transitions  (1-3).  The  most  recent  comprehensive  studies  of  these  bands 
were  reported  by  Wicland  (4-6)  over  20  years  ago.  Wieland’s  HgBr  and  Hgl  sources 
contained  HgXj  in  the  natural  isotopic  mixture.  In  the  case  of  HgCI  he  used  isotop¬ 
ically  enriched  CI2,  but  still  natural  Hg.  To  refine  the  spectroscopic  characterization 
of  these  systems,  we  have  undertaken  a  study  of  the  B  X  transitions,  using  iso¬ 
topically  pure  mercury  and  halogen  in  our  sources.  Such  single-isotope  sources  have 
been  used  to  great  advantage  in  previous  works  from  this  laboratory  (7-9). 

We  have  recently  reported  preliminary  results  of  our  work  on  HgBr  (/(?),  and  HgCI 
and  Hgl  (11).  In  this  paper  we  present  our  detailed  analysis  of  the  B  X  transition 
in  Hgl — the  band  with  peak  intensity  near  4450  A.  Our  interpretation  of  the  spectrum 
agrees  qualitatively  with  Wieland's  work  (6),  except  at  the  long  wavelength  end. 
However,  our  least-squares  analysis  indicates  that  the  u*  numbering  suggested  by 
Wieland  is  one  unit  too  high.  Our  computational  analysis  includes  band-profile  sim¬ 
ulations  and  Franck-Condon  calculations  to  corroborate  the  assignments  of  many 
of  the  features  in  the  spectrum  and  to  deduce  the  relative  configuration  of  the  B  and 
X  potential  curves. 

EXPERIMENTAL  SECTION 

It  recording  the  spectra,  we  used  procedures  similar  to  those  described  previously 
(7-9).  The  sources  consisted  of  3-mm  o.d.  pyrex  tubes  about  10  cm  long,  which  were 
charged  initially  with  ^Hg'^^I:  (or  ’•’’’Hg'”!:)  and  ~2(X)  Torr  Ar,  then  sealed  off 
with  a  torch.  The  ^Hg  (95.7%  isotopic  purity.  Oak  Ridge)  was  obtained  in  the  form 
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by  our  analysis.  ■  -  2  sufficed  for  both  tbe  B  state  and  tbe  D  state.  constructing  and  solTlng  the  equations  for  the  polynonlal  and  NDE  fits  hare 


thu  did  iht  polyooalals  (r  .  Fig.  2).  tbca  lead  powers  were  dropped  froa  tbe 


laposea.  raere  was  aegiigioie  lacrease  in  rariance.  or  course  tnese  alxed  fits  failed  to  give  satlsfactorj  approach  to  dissociation,  because  they  "turned  on" 
have  sore  adjustable  paraaeters  than  do  the  single  representation  fits,  so  they  too  close  to  v^.  Here  we  define  "satisfactory*  behavior  as  qualitative 
are  less  efficient.  However,  they  have  soae  advantages  which  are  worth  the  agreeaent  with  the  graphical  procedure  which  served  as  a  progenitor  of  the 
price  of  the  additional  paraaeters.  as  discussed  below.  direct  fitting  aethods  [16.19].  The  graphical  aethod  consists  of  extrapolating 

the  highest  observed  levels  to  dissociation  using 


reliable  for  other  Isotopoaers  of  HgBr.  through  use  of  the  staadard  Isotopic 
relations  [13.171.  Ve  are  currently  trying  to  «siend  the  B*X  asslgnnents  to 
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Appendix  7 


•Br5T’  SPECTROSCOPIC  CONSTANTS  FOR  HgBr  ^ROM  DIRECT  FITS  OF  MULTIPLE  BAND  SYSTEMS  TO 
POLYNOMIALS  AND  NEAR-DISSOCIATION  EXPANSIONS 

J  GAII.  ASHMORF  AND  JOEL  TELLINGHUISEN 

The  B-X  (4200-5100  A).  C-X  (2700-2950  A),  and  D-X  (2480-2700  A)  transitions  OT  Hpr  have 

200  79  200  81 

been  photographed  and  analyzed  for  Isotoplcally  pure  Hg  Br  and  Hg  Br.  The  analyses 
yield  laproved  vibrational  constants  for  all  four  states  and  rotational  constants  for  the  B 
and  X  states.  Optlaal  spectroscopic  paraneters  arc  obtained  for  all  four  states  froa  direct. 
siBUlianeous  fits  of  all  three  transitions  to  the  standard  polynoaials  in  (v*l/2)  and  to 

1.2 

near-dissociation  expansions. 

In  addition  to  the  above-mentioned  systems,  we  have  recorded  and  analyzed  by  computer 
siBOlaiion  the  B-A  transition  (5500-8000  X).  Efforts  are  currently  underway  to  (1)  measure 
collisional  line  broadening  in  the  B-X  system  using  a  Fabry-Perot  interferometer,  and  (2) 
determine  the  R-dependence  of  the  B-X  transition  strength  function  from  analysis  of  relative 
intensity  data. 


*  Work  supported  by  the  Office  of  Naval  Research. 

*R.  J.  Le  Roy  and  W-H.  Lam.  Chem.  Phys.  Lett.  JX-  544  (1980). 
"J.  Tellinghuisen.  J.  Chem.  Phys.  (in  press). 
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interfacing  a  MICRODENSITOMETER  TO  A  MICROCOMPUTER* 
r  CAR!  YSLE  SALTER  AND  JOEL  TELLINGHUISEN 

In  methods  of  photographic  spectroscopy  there  is  a  need  for  precision  measurement  of  a 
large  amount  of  experimental  data  —  the  positions  and  intensities  of  rotational  lines 
vibrational  band  heads,  and  calibration  lines  on  the  photographic  plate.  Microdensltometers 
and  optical  comparators  permit  one  to  measure  the  positions  of  sharp  features /^th  a 
precision  of  1-2  pa:  however,  the  procedure  of  measuring,  recording,  and  logging  the  feta  for 
further  computer  processing  can  be  very  tedious  and  time  consuming,  if  done  manually  To 
expedite  this  aspect  of  our  work,  we  have  designed  and  built  a  cheap  (-$2500.  microcomputer 
included)  control  interface,  by  means  of  which  a  TRS-80  Model  III  microcomputer  controls  the 
motion  of  the  plate  on  a  mlcrodensitometer  and  logs  the  optical  density  in  digital  form  In 
this  paper  we  di.scuss  various  aspects  of  the  interfacing  task,  including  hardware  and 
software  for  stepping  motor  control,  analog-to-dlgltal  conversion,  and  extraction  of  line 
positions  from  the  recorded  data. 


*  Work  supported  by  the  Office  of  Naval  Research. 

ill  AaVtigrs :  Department  of  Chemistry, 

Nashville.  TN  37235. 
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